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Abstract  This  paper  addresses  the  problem  of  beamforming  a  uniform  linear  array 
when  some  of  the  receive  elements  are  not  operational.  Bad  sensors  are  often 
handled  by  either  zeroing  or  interpolating  the  faulty  elements  prior  to  conventional 
beamforming.  While  zeroing  faulty  elements  prior  to  conventional  beamforming  is 
the  simplest  approach,  it  often  results  in  undesirably  high  beamformer  sidelobes. 
Alternatively,  minimum  mean-square  error  (MMSE)  interpolation  of  the  missing  data 
is  not  explicitly  aimed  at  minimizing  post-interpolation  leakage  of  strong  interference 
components  into  otherwise  quiet  directions.  While  true  minimum  variance  (MV) 
adaptive  beamforming  is  the  optimal  solution  given  long  observation  times,  for  large 
arrays  in  highly  dynamic  environments,  severely  limited  snapshot  support  poses  a 
difficult  trade-off  between  desired  interference  suppression  versus  unwanted  signal 
cancellation. 

In  this  work,  we  propose  an  alternative  approach  for  conventional  beamforming  with 
faulty  sensors  based  on  adaptively  synthesizing  complete  array  data  snapshots  which 
minimize  post-interpolation  leakage  into  quiet  directions  subject  to  a  constraint  that 
the  solution  is  sufficiently  close  to  the  data  measured  at  the  working  elements. 
After  reconstruction  of  the  complete  array  data  snapshots,  computationally  efficient 
conventional  beamforming  can  be  performed  to  both  estimate  the  noise  field 
directionality  as  well  as  produce  time-series  output  for  further  temporal  analysis. 
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Beamforming  Sensor  Arrays  with  Faulty  Sensors 


OBJECTIVE:  To  mitigate  sidelobe  degradation  due  to  sensor  outages  in  large  sonar  arrays. 

BACKGROUND: 

•  Target  detection  in  interference-dominated  littoral  regions  depends  critically  on  achieving  low  beamformer 
sidelobe  levels.  Sensor  failures  can  seriously  degrade  the  array  pattern  of  conventional  beamformers. 

•  While  zeroing  faulty  elements  prior  to  conventional  beamforming  is  a  common  remedy,  it  often  results  in 
undesirably  high  beamformer  sidelobes. 

•  Adaptive  beamforming  using  only  the  good  sensors  provides  asymptotically  optimal  array  gain,  however 
highly  dynamic  littoral  environments  make  achieving  this  asymptotic  performance  a  challenging  problem. 

•  We  investigate  adaptive  channel  interpolation  and  compensation  techniques  which  mitigate  faulty  sensors 
on  a  snapshot-by-snapshot  basis  and  are  thus  applicable  to  large  arrays  in  dynamic  environments. 

•  Minimum  mean-square  error  (MMSE)  interpolation  of  bad  elements  is  an  obvious  approach  but  doesn’t 
explicitly  minimize  leakage  of  strong  interference  components  into  otherwise  quiet  directions. 

•  Spatial  autoregressive  (AR)  methods  have  been  previously  proposed  for  array  interpolation  by  Swingler 
and  Walker  (1988)  but  are  not  designed  to  handle  spatially  inhomogeneous  multipath  interferers. 

•  In  this  work,  we  consider  a  linear  inverse  theoretic  framework  for  mitigating  the  effects  of  faulty  sensors. 
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Motivation  for  Mitigating  Faulty  Array  Elements 


•  Bearing-time  record  (BTR)  from  real  30-element  passive  sonar  towed-array 
data.  Full-array  BTR  (left).  Degraded  with  6  zeroed  channels  (right). 

•  Note  degradation  in  sidelobe  response  complicates  weak  target  detection. 


sin[B) 
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Array  Data  Modeling  with  Missing  Sensors 


•  Model  the  output  from  an  7V-element  uniform  linear  array  as  zero-mean  complex 
snapshots  xA.  with  covariance  matrix  consisting  of  signal,  interference,  and  noise: 

Rx=a;d(ds)d(0sy+AVA++a-l 

where  d (6 s )  and  o s  are  the  signal  wavefront  and  power,  APA+  defines  the  unknown 
interference  covariance,  and  <Jn\  represents  the  white  noise  component. 

•  Zeroing  missing  sensors  corresponds  to  using  a  covariance  matrix  with  correspondingly 
zeroed  columns/rows  e.g. 
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•  Classical  interpolation  methods  can  be  thought  of  as  “filling  in”  missing  covariance 
matrix  elements  using  Toeplitz  approximation  and  then  designing  a  linear  filter  for 
MMSE  interpolation  on  a  snapshot-by-snapshot  basis. 
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Covariance  Matrix  Estimation  with  One  Snapshot 


•  Covariance  matrix  estimation  with  as  little  as  a  single  snapshot  of  data  typically 
requires  imposing  a  Toeplitz  structure  (i.e.  spatially  stationary  field). 

•  A  straightforward  Toeplitz  covariance  matrix  estimate  can  be  computed  by  using  the 
weighted  projected  covariance  estimator  proposed  by  Barton  and  Smith  (1997)  which 
uses  biased  average  diagonals  of  the  sample  covariance  matrix. 

•  For  a  fully  populated  array,  this  method  is  equivalent  to  computing  a  Toeplitz 
correlation  matrix  with  the  biased  autocorrelation  sequence  of  the  array  data. 


n  =  k 


•  Toeplitz  approximation  models  the  interference  subspace  as  a  linear  combination  of 
uncorrelated  plane  wave  returns  corresponding  to  interference  directions: 


and  P  and  Amin  are  positive  diagonal  matrices. 

•  More  sophisticated  methods  for  structured  covariance  matrix  estimation  have  been 
developed  by  numerous  researchers  (e.g.  Jansson  and  Ottersten,  2000,  Robey  and 
Fuhrmann,  1991,  Stoica  and  Li,  1999,  Cadzow,  1988,  Williams  and  Johnson,  1988,  ) 
but  typically  involve  greater  computational  complexity. 
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Minimum  MSE  Interpolation  of  Faulty  Sensors 


A  natural  solution  for  handling  faulty  channels  is  to  interpolate  them  using  a  linear 
combination  of  the  good  sensor  outputs. 

The  linear  filter,  T/? ,  which  minimizes  mean-square  error  (MSE),  i.e. 

2 


min 


Xb-TbXg 


is  given  by  Tb  =  Rb  R 


-1 

bg^'gg 


where  x  =  L+x  are  the  N-m  good  channels  and  xb  —  K+x  are  the  m  bad  channels. 


•  A  single  snapshot  estimate  of  the  cross-correlation  between  faulty  and  good  channels  is 
obtained  from  the  Toeplitz  covariance  matrix  estimate  of  the  full  array  data. 

•  The  Minimum  MSE  (MMSE)  interpolated  full  array  data  snapshot  can  be  written  as  a 
linear  transformation  matrix  Tm  operating  on  the  original  data  snapshots: 

qra  =  TmX  where  Tm  =  (L  +  KTb)LH 


•  Minimizing  channel-wise  MSE  does  not  explicitly  minimize  sidelobe  leakage  of 
interferers  into  quiet  directions  and  thus  only  indirectly  mitigates  impact  of  bad  channels. 

•  Alternative  ALPINEX  method  of  uses  AR  modeling  of  interferers  and  averaged  forward 
and  backward  linear  prediction  of  faulty  channels  (Swingler  and  Walker,  1998). 
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Adaptive  Channel  Compensation  (ACC) 


Idea:  Adaptively  reconstruct  full-array  snapshot,  q,  such  that  interference  leakage  into 
nominally  quiet  directions  is  minimized. 


The  proposed  adaptive  channel  compensation  (ACC)  method  minimizes  interference 
leakage  into  quiet  plane-wave  directions  subject  to  the  constraint  that  the  data  distortion 
at  the  good  elements  is  within  some  tolerance,  £  i.e. 


min  (q  <M>  q)  subject  to  the  constraint 
q 


LHq-xs 


<£ 


Solution:  q  =  (<M>+  +  //  LL+ ) 1  ju  Lxg 

where  the  columns  of  O  are  the  sub-dominant  eigenvectors 

/\ 

of  the  Toeplitz-approximated  Rx. 


•  Choice  of  £  (which  is  inversely  related  to  A)  determines  the  trade-off  between  signal 
wavefront  distortion  and  mitigation  of  interference  leakage. 
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ACC  as  the  Solution  to  a  Linear  Inverse  Problem 


•  Adaptive  channel  compensation  can  be  interpreted  as  inversion  for  the  complete  data 
using  the  linear  operator  T  operating  on  the  corrupted  data  snapshot: 

T  =  (4><t»+//  LL)'//LL 

•  Decomposing  T  in  terms  of  the  dominant  ^max^^axand  sub-dominant  00+  subspaces 
of  Rx  gives: 

T  =  +  (-^-)<&®+  )(I  -  K(K+«DOK) 1  KM*  )  LL+ 

1  +  n 

where  K+  and  L+  select  the  faulty  and  good  elements  and  LLH  +  KKH  =  I 

•  Effect  of  ACC  is  to  make  interferers  look  more  like  plane-waves  across  the  full-array 
which  can  then  be  effectively  nulled  by  conventional  beamformer  shading. 

•  Signals  in  the  subdominant  subspace  are  not  interpolated  but  suppressed  by  a  factor  of 

JLl/(\  +  JLl) 

•  Choice  of  jd  =  1 0  empirically  found  to  provide  good  compensation  of  faulty  channels 
with  negligible  signal  distortion. 
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Comparing  ACC  vs.  Adaptive  Beamforming 


•  Owsley  (1984)  previously  studied  asymptotic  array  gain  (AG)  improvement  achieved  by 
adaptive  beamforming  versus  conventional  processing  in  terms  of  Rg  =  AG0ptimum  j AGCB 


Using  this  metric  to  compare  conventional  beamforming  with  ACC  vs.  adaptive 
beamforming  (ABF)  gives: 


SINR0 
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where  L, 


d(Y?  )H  Td(  ^ ) 
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For  sidelobe  interference,  LtU  1  ,  and  thus 
ACC  only  slightly  under-performs 
asymptotically  optimum  beamformer  for 
moderate  INR’s. 

ACC  expected  to  have  the  advantage  in  non- 
asymptotic  scenarios  where  ABF  nulling 
degraded  by  poor  eigenvector  estimation. 


Array  gain  improvement  ratio  for  optimum  MVDR  beamformer 
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Example  Beamformer  Array  Pattern  Degradation 


•  Simulation  example 
of  the  array  patterns 
(right)  from  a 
Hamming-shaded 
30-element  array 
data  with  and 
without  6  missing 
elements  (below). 
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•  Note  severe  sidelobe 
degradation  due  to 
zeroed  elements. 
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Simulated  Single  Snapshot  Field  Directionality  Map 


•  Comparison  of  clairvoyant  full  array,  zeroed  conventional,  ALPINEX,  and  MMSE 
methods  versus  adaptive  channel  compensation  (ACC)  for  a  30  element  array 
with  a  0  dB  target,  20  dB  correlated  interference  at  sin(0)  =  0  and  -0.2  and  10 
dB  uncorrelated  plane-wave  interference  at  sin(0)  =  .47 


\ 


ACC  with  Hamming  beamforming  clearly  unmasks  the  target  and  achieves  lowest 
sidelobe  levels  for  this  scenario  of  6  faulty  sensors. 


-  Full  Array  Conventional 

-  Array  w/  faulty  elements 

-  ACC 

-  ALPINEX 

-  MMSE 
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ACC  vs.  ABF  with  Limited  Sample  Support 


SINR  vs.  snapshot  support  simulation  for  30  sensor  array  with  6  missing  sensors  in  simulated 
multipath  environment  and  target  SNR  =  -10  dB. 

ACC  comparable  with  dominant-mode  rejection  DMR-MVDR  (Cox,  Owsley)  and  superior  to 
MMSE  and  SMI-MVDR  for  small  sample  support. 

Unlike  all  the  ABF  methods,  ACC  useful  down  just  one  snapshot. 
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SINR  as  a  function  of  number  of  snapshots 
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—  Clairvoyant  Full  Array 

—  Array  w/  faulty  sensors 

—  ACC 

—  MMSE 

—  MVDRw/DMR 
MVDR 

MVDR  full  array 
--  DMR  full  array 
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ACC  Performance  vs.  Conventional  Interpolation 


•  Single  snapshot  ROC  for  -4  dB  target  (30  sensor  array  with  6  missing)  in  simulated  multipath 
environment  (left).  PD  vs.  SNR  @  PFA=0.04  (right). 


I 


•  ACC  slightly  outperforms  MMSE  at  all  SNR’s. 


•  ALPINEX  (Swinger  and  Walker)  degrades  significantly  in  correlated  multipath. 


Receiver  Operating  Characteristic  (SNR  =  4  dB) 


Probability  of  detection  as  a  function  of  SNR  for  PF  =  0.04 
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Simulated  BTR  for  ACC  vs.  Clairvoyant  Beamformer 


Bearing  Time  Record  -  ACC 
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•  Bearing-time  records  (BTR’s) 
counterclockwise  from  above:  Single¬ 
snapshot  ACC,  Clairvoyant  CBF,  and 
CBF  with  zeroed  elements.  Target  at 

sin(0)  =  0.27  @  0  dB  SNR. 

•  ACC  clearly  unmasks  target  in  the 
presence  of  resolvable  interference. 
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Simulated  ACC  vs.  Interpolated  Beamformers 


Bearing  Time  Record  -  ALPINEX 
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•  Bearing-time  records  (BTR’s) 
counterclockwise  from  above:  Single¬ 
snapshot  ACC,  ALPINEX,  and  MMSE 
CBF  with  zeroed  elements.  Target  at 

sinlT?)  =  0.27  @  0  dB  SNR. 

•  ACC  sidelobe  suppression  better  than 
ALPINEX  with  fewer  “ghost”  tracks 
apparent  in  MMSE  solution. 
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Real  Data  BTR  for  ACC  vs.  Clairvoyant  Beamformer 


Bearing  Time  Record  -  Full  Array 
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Bearing  Time  Record  -  ACC 
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•  Bearing-time  records  (BTR’s) 

counterclockwise  from  above  for  real 
30-element  towed  array  data:  Single¬ 
snapshot  ACC,  Clairvoyant  CBF,  and 
CBF  with  zeroed  elements. 
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Lack  of  ground-truth  limits  the  ability  to 
draw  definitive  conclusions. 
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Real  Data  ACC  vs.  Interpolated  Beamformers 


Bearing  Time  Record  -  ALPINEX 
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Bearing  Time  Record  -  MMSE 
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Bearing  Time  Record  -  ACC 
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•  Bearing-time  records  (BTR’s) 
counterclockwise  from  above 
for  real  30-element  towed  array 
data:  Single-snapshot  ACC, 
ALPINEX,  and  MMSE  CBF. 


ACC  comparable  to  MMSE  in 
this  scenario. 
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Summary  and  Further  Work 


•  With  increasing  array  apertures  and  channel  count,  mitigating  the  sidelobe  impact  of  bad 
elements  may  become  a  critical  problem. 

•  ABF  snapshot  support  problem  exacerbated  with  large  arrays  in  dynamic  littoral 
environments. 

•  MMSE  interpolation  methods  not  explicitly  aimed  at  reducing  leakage  of  strong 
interference  into  quiet  look  directions. 

•  ACC  proposed  as  a  means  of  compensating  the  entire  array  so  that  strong  interferers 
made  to  look  more  like  uncorrelated  plane-waves  which  can  then  be  suppressed  by 
conventional  beamformer  shading. 

•  ACC  outperforms  MMSE  and  ALPINEX  in  terms  of  output  SINR  in  low-snapshot 
support  correlated  multipath  environments. 

•  Further  work  required  to  study  performance  in  the  presence  of  sensor  calibration  errors. 

•  ACC  recently  applied  to  372  element  OTH  radar  data  with  success  in  the  presence  of 
faulty  elements. 

•  Possible  applications  of  ACC  to  sonar  data  might  include  TB-29  or  large  bottom- 
mounted  horizontal  arrays. 
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